PERSISTENCE PROBABILITIES & EXPONENTS 

FRANK AURZADA AND THOMAS SIMON 

^^ ' Abstract. This article deals with the asymptotic behaviour as i ^- +00 of the survival 

^— V , function F[T > t], where T is the first passage time above a non negative level of a random 

rv^ ' process starting from zero. In many cases of physical significance, the behaviour is of the 

J, , type F[T > t] ^ ^-^+0(1) for a known or unknown positive parameter 6 which is called a 

C^ ' persistence exponent. The problem is well understood for random walks or Levy processes 

but becomes more difficult for integrals of such processes, which are more related to physics. 

We survey recent results and open problems in this field. 
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1. Introduction 

Let {Xt, t > 0} be a real stochastic process in discrete or continuous time, starting 
from zero. The analysis of the first passage time Tx = inf{t > 0, Xt > x} above a non- 
negative level X is a classical issue in probability. In this paper we will be concerned with 
the asymptotic behaviour as t — ?• +00 of the survival function F[Tx > t] for a class of 
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processes related to random walks and Levy processes. This problem has attracted some 
interest in the recent literature under the denomination persistence probability. In the 
self-similar framework, it is also related to the lower tail probability problem which is the 
asymptotic study as e — )■ of the quantity F[St < e] where St = sup{Xs, s < t} is the 
one-sided supremum process. In many situations of interest it turns out that the behaviour 
is polynomial: one has 



;i.i) p[T. >t] = f 



-e+o{i) 



for a non- negative parameter 6 called a persistence exponent, which usually does not depend 
on x > and often belongs to [0, 1]. The study of asymptotic behaviours like (11.11) has gained 
some attraction over the last years in the physical literature as well, where the parameter 6 
is often called a survival exponent. Notice that estimates of the type f ll.ip also appear in 
reliability theory, a subject that we shall here not discuss, where T^ is viewed as a certain 
failure time whose typical upper tails are Pareto-like - see [5^ . 

If X is a random walk or a Levy process, studying the law of T^. is a special part of 
fluctuation theory and the persistence probabilities are then well understood. For example 
if Xi is attracted to a stable law, classical fluctuation identities entail that the persistence 
exponent is the positivity parameter of the latter. There are many accounts on fluctuation 
theory and we refer in particular to [9j for random walks and to [20j for Levy processes. In 
the first part of this article we focus on the results of this theory dealing with the asymptotic 
behaviour of P[Tj,. > t], and we try to be as exhaustive as possible. In general, the behaviour 
is the same in discrete and continuous time, except of course for a; = where the problem 
becomes different for Levy processes. We notice that even though a posteriori the resulting 
exponents turn out to be the same, there is no simple approximation argument that would 
yield this a priori, and for this reason we have to consider the discrete and continuous time 
situations separately. The results, all classical, are presented here in order to give some 
insight into more complex situations where X, though constructed upon a Levy process or 
a random walk, is not a Markov process anymore. 

Some of these more complex situations, which are called non-trivial in the physical liter- 
ature, are the matter of the second part of this article. We first consider integrated random 
walks and integrated Levy processes. It turns out that for such simple constructions, the 
computation of the persistence exponent is not quite easy in general. We display recent and 
less recent results where the persistence probabilities are estimated with various degrees of 
precision, and which all have the common feature that the persistence exponent of the under- 
lying process is twice the persistence exponent of its integral. We believe that this is a kind 
of universal rule and several conjectures are stated in this direction, the most general being 
Conjecture 3 and Conjecture 5. A curious fact is that the situation with no positive jumps, 
where classical fluctuation theory becomes much simpler, appears to be more difficult than 
the dual situation with no negative jumps, where some exponents of integrated processes 
can be computed. The reason for that seems to be only methodological, and to cast out the 
curse of spectral positivity in the study of persistence probabilities of integrated processes 
probably requires other, say deeper, tools. 
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We then consider fractionally integrated processes, where the situation is more difficult. 
The case when the underlying Levy process is Brownian motion yields self-similar Gaussian 
processes which can be changed into a stationary one by the Lamperti transformation. The 
sought after persistence exponent is then directly related to an estimate on the probability of 
non-zero crossings for these Lamperti transforms, a problem known to be hard in spite of the 
Ito-Rice formula which gives some information at the expectation level in the smooth case. 
We present some universality and monotonicity results, and also some partially heuristic 
comparisons with the fractional Brownian motion case, whose persistence exponent can be 
computed explicitly. We also display other explicit computations of persistence exponents 
for related processes such as weighted sums or iterated processes. 

The third and last part of this article deals with some applications of the persistence 
probabilities in mathematical physics. We ffist deal with Lagrangian regular points of the 
inviscid Burgers equation with random self-similar data. The link between the Hausdorff 
dimension of such points and the persistence exponent of integrated processes dates back to 
the original paper [52] . We recall here that the problem is still open in the fractional Brownian 
motion case and state a plausible conjecture when the initial data is a two-sided stable Levy 
process. Second, we consider the zero-crossings of a peculiar Gaussian stationary process 
which is related to the positivity of Kac polynomials with large even degree. This connection, 
which was discovered in [T7], has further ramifications with the persistence exponents of 
integrals of Brownian motion with higher order and of a certain diffusion equation with 
white noise initial conditions in the plane, and we make a brief account on the subject. We 
last consider three different interacting statistical systems whose analysis hinges significantly 
upon the persistence probabilities of integrated processes: wetting models with Laplacian 
interaction, fluctuating interfaces with Langevin dynamics, and sticky particles on the line 
with Poissonian initial conditions. 

Some open problems stated in the present paper are believed to be challenging and we 
think that they could catch the attention of some colleagues. We finally point out that we 
have not exhausted here all implications of persistence in physics and that the persistence 
exponent of many other models remains unknown - see fTEl [55] and the references therein. 

2. Classical results 

2.1. Random walks. Let {X„, n > 1} be a sequence of i.i.d. real random variables with 
common distribution fi, and Sn = Xi + . . . + X„ be the associated random walk. Consider 

Tj. = inf {n > 1, S'„ > x} 

the first-passage time above a; > 0. Recall the following basic rule, to be found in e.g. Chapter 
XILl-2 of [221 : 

(2.1) Tj. is a.s. finite for every x > <^ Tq is a.s. finite, 

and that the latter is also equivalent to the fact that Sn does not drift to — oo. In this case 
one has P[Tj. > n] — )■ as n — > -|-cx) for every x > 0, and the difficulty to estimate the 
speed of convergence comes from the fact that the event {T^ > n} = {Si < x, . . . , Sn < x} 
depends on n correlated random variables. 
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If n is concentrated on M"*", then P[Ta. > n] = F[Sn < x] and the problem becomes one- 
dimensionaL The straightforward inequahty 

shows that F[Tx > n] tends to zero at least exponentially fast (unless /i is degenerate at zero). 
Notice that the rate might also be superexponential and depend on x at the logarithmic scale: 
if Xi has a positive strictly a— stable law for instance, de Bruijn's Tauberian theorem - see 
Theorem 4.12.9 in [10] - leads to 



-logP[T^>n] = -logP[5„<a;] = -\ogF[Si < xn-^^''] ~ k^x-"/^^-")^^/^^-") 

for every a; > 0, with some explicit /€„ > 0. Of course, one has Tq = 1 a.s. whenever /i does 
not charge 0. If /i{0} > 0, Jain and Pruitt's general uniform results on renewal sequences 
express the asymptotic behavior of F[Sn < x] in terms of quantities related to fi, at the 
logarithmic scale - see Theorem 2.1 in |30] - and also at the exact scale, under some conser- 
vativeness assumption - see Theorem 4.1 in [30] . 



If n is not concentrated on M"*", it is easy to see that the asymptotics of P[T^ > n] will not 
drastically depend on x: choosing e > such that /i^ = /i(— oo, —e) > 0, the Markov property 
entails P[T,. > n] > P[To > n] > F[T, > n- l]P[Xi < -e] > I2,F[T, > n] > /i^P[T, > n] as 
soon as ke > x, so that 

P[T, > n] X P[To > n] 

for every x > 0. 

The following formula computes the generating function of {P[To > n], n > 0} in terms 
of that of the sequence {n^^F[Sn < 0], n > 1}, and is true for any random walk. 

Sparre Andersen's formula. For every 2; g] — 1, 1[ one has 



(2.2) ^2"P[To>n] = exp 



n>0 






n>l 



This result is obtained after simple rearrangements from Theorem XII. 7.1 in [22j, whose 
proof has a combinatorial character and depends heavily on the independence and stationar- 
ity of the increments of the random walk. A simpler method relying on elementary Fourier 
analysis - see Chapter 3.7 in [21] - yields the more general 

Spitzer's formula. For every 5; g] — 1, 1[ and A > 0, one has 






-A5+1 



n 

.n>l 



(2.3) E2"E[e-^^"] = exp 

n>0 

with the notation M„ = max(0. Si, ... , S'„) and S*^ = max(0, Sn). 

Indeed, one obtains (12. 2p as a consequence of (12. 3p in letting A — )■ +00. It is beyond the 
peculiar scope of the present paper to discuss the full strength and the various generalisations 
of Spitzer's formula such as Baxter-Spitzer's formula or the Wiener-Hopf factorization, and 
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we refer e.g. to [91|201|22] for more on this subject. Let us simply remark that Sparre Ander- 
sen's formula entails easily - see e.g. Theorem XII. 7. 2 in [22] - the following characterization 
of (J2l]): 

Sn drifts to -oo <^ ^^ - P[S'„ > 0] < +oo. 

n>l 

From now on, we will suppose that Sn does not drift to — oo and that /i is not concentrated 
on M+. This entails that P[S'„ > 0] > and P[S'„ < 0] > for every n > 1. 

A remarkable consequence of (12.21) is that when W[Sn < 0] = p G (0, 1) for every n > 1, 
one obtains an explicit formula for P[To > n] depending only on p and n : one has 

^ ^ " ^ (l-z)p ^ n!r p 

n>0 ^ ' n>0 ^^' 

SO that 

(2.4) F|r„ > n\ = ^-^±^ ~ 'f-. 

^!r(p) r(p) 

For example, symmetric random walks such that P[5'„ = 0] = for every n > 1 (this latter 
property is true when [i is non atomic, for instance) all enjoy the property that 

r(n + i/2) 1 



P[To > n\ 



y^n\ a/v™ 



Recall in passing that the estimate is slightly different for simple random walks with n{l} 
p{—l} = 1/2, since then F[S2n = 0] 7^ 0. The classical computation 



and the Tauberian theorem for monotonic sequences entail P[To > n] ~ a/2/7™. 

Another remarkable consequence of (12.31) is the exact computation of the persistence ex- 
ponent whenever {Sn,n > 1} fulfils the so-called Spitzer's condition 

(2.5) lim F[Sn < 0] = pe (0,1). 
A theorem of Rogozin shows indeed that (12.51) entails 

(2.6) P[r. > n] ~ ^^!QM = ^p-i+o(i) 

for every x > with l(n) some slowly varying sequence and Cx some explicit positive constant. 
Besides, as explained in Theorem 8.9.12 of ^Q\, this asymptotic behavior is actually equiva- 
lent to (12. 5p . Recall - see Chapter 7 in |20] - that (12. 5p is also equivalent to P[S'„ < 0] — )• p 
in Cesaro mean. When F[Sn < 0] — )■ various behaviors are possible, contrary to the above. 
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For example if nPF[Sn < 0] is slowly varying for some p G (0, 1), Theorem 8.9.14 in ^0] and 
Theorem XVII.5.1 in [22] yield 

(2.7) P[To >n] ^ l{n)n-^^+P^ 

with l{n) some slowly varying sequence. Other various behaviours also appear when p has 
positive expectation and we refer to [I9] for precise results. When P[S'„ < 0] — ;■ 1, there does 
not seem to be any universal behavior for W[Tq > n], although one might expect that the 
persistence exponent is always zero. 

We conclude this paragraph with the so-called upward skip- free or right-continuous random 
walks on Z viz. such that Supp fi C {1, 0, —1, —2, —3, . . .}, for which our survival analysis 
does not require the use of (12.21) . Indeed, the distribution of T^ is then given by Kemperman's 
formula [33j which reads: 

(2.8) F\Tk = n] = —^ FlSn = k + 1] for every n> k. 

n 

In this particular case, we are reduced to the asymptotical behavior of the one-dimensional 
probability P[S'„ = 1]. Notice - see |3l] - that there is a universal upper bound P[5'„ = 
1] < cn-^/2 leading to P[Tfc > n] < 2c{k + l)n-i/^ but the exact behavior of F[Sn = 1] 
depends on /i. For example, if fi is in the respective domain of normal attraction of some 
strictly a— stable law - with a G (1,2] since otherwise S'„ would drift to — oo by the law 
of large numbers, Gnedenko's local limit theorem - see e.g. Theorem 8.4.1 in [lOJ - yields 
F[Sn = 1] ~ Can~^f°'l{n) for some explicit Cq,, so that by (12. 8p 

(2.9) P[Tfc > n] ~ acc,{k + l)l{n)n~^''' . 

Recall - see Proposition 8.9.16 in [10] - that p is in the domain of attraction of some strictly 
a— stable law with a G (1, 2] if and only if (12. 5 p holds with p = 1 — 1/a, so that (12. 9 p is 
actually a consequence of (12. 6p . On the other hand, the slowly varying term l{n) can be 
removed if Sn is in the respective domain of normal attraction (viz. when Sn/n^^"' converges 
in law to some non degenerate limit - see e.g. the concluding remark of Chapter XVII. 5 in 
and this degree of precision is not given by Rogozin's theorem. 



2.2. Levy processes. Let {Zt-, t > 0} be a non-degenerate real Levy process starting from 
and 

T, = inf{t > 0, Zt > x} 
be its first passage time above a; > 0. Consider {Z„, n > 1} the associated random walk. 
The inequality 

(2.10) p[r, >t\ < p[r, > \t\] 

with the notation T^ = inf {n > 1, Zn > x} , yields a rough upper bound for P[Tj; > t] which 
can be made more precise as a function of t and x in applying the results of the previous 
paragraph. This upper bound however does not yield enough information in general if x = 0. 
For example Rogozin's criterion - see Proposition VI. 11 in [S] - shows that Tq = a.s. when 
Z has unbounded variations, whereas the function t i— )■ F[Zi < 0, . . . , Z^t] < 0] has a positive 
limit at +oo if Z drifts to — oo. Recall also that in the bounded variation case, the regularity 
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of the half-line for Z is characterized in terms of the Levy measure - see Theorem 22 in 
and again P[Zi < 0, . . . , Z^j] < 0] might have a positive limit when t — )■ +oo eventhough 
To = a.s. 

If X > however, it turns out that the two quantities in fl2.10p are often comparable. First 
of all, for every x > one has 

(2.11) P[T^ = +oo] > ^ Zi ^ -oo a.s. ^ / - P[Zi > 0]dt < +oo 

(see e.g. Proposition 4.6 in [20]), and a simple analysis shows that the latter condition is 
equivalent to 



Zn -^ -OO a.s. ^ V - Wn > 0] < 
< ^ r) 



n>l 



OO, 

n 



so that with the above notation F[Tx > t] — )■ if and only if F[Tx > t] — )■ as t — ?■ +oo. In 
the following, we will see that in many examples one has 

P[T, > t] X F[f, > [t]]. 

However, it does not seem easy to prove this estimate a priori, which probably does not hold 
in full generality. 

We will suppose henceforth that (12. lip does not hold. Set /i for the law of Zi, Mf = 
sup {Zs,s < t} for the running maximum of Z, and recall that 

P[T, >t] = F[Mt < x]. 

If fi is concentrated on M+ then Z is a.s. increasing and one is reduced to the random walk 
case because 

P[^[t]+i < a;] < P[T, > t] = F[Zt < x] = F[Mt < x] < ¥[Z[t] < x]. 

If /i is not concentrated on M+ then an argument analogous to that of the previous paragraph 
shows that 

P[T, > t] X ¥[Ty > t] 

as t — )■ +00 for every x,y > 0, as in the discrete time framework. The classical approach 
to obtain more information on F[Tx > t] relies on a particular case of the first fluctuation 
identity - see e.g. Theorem VI. 5 in |5], which is an analogue of Spitzer's formula in continuous 
time: 

Baxter-Donsker's formula. For every A > and q > 0, one has 

"OO -gt 



(2.12) E[e"^'"^'] = exp 






with the notation Z^ = max(0, Z^) and Tq ~ Exp (g) an independent random time. 



-oo 
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An important consequence of this formula is a theorem of Rogozin - see e.g. Theorem 
VI. 5 in [5] - which shows that if the following Spitzer's condition 

(2.13) - / F[Zs>0]ds -^ pG (0,1), t^ 
holds, then for any x > 

(2.14) p[r, > t] ~ c,/(t)t-^ 

with Cx > and l(t) some slowly varying function not depending on x. Besides, one can 
show that (12.131) and (I2.14p are actually equivalent - see again Theorem VI. 5 in [5]. Re- 
call - see Chapter 7 in j20] - that (I2.13P is also equivalent to (12. 5 p for the random walk 
{—Zn,n > 1}. The estimate (I2.14p can be refined for strictly a— stable processes, which all 
enjoy the property that 

(2.15) F[Zt > 0] = p G (0, 1) for every t > 0. 

An asymptotic analysis of the so-called Darling integral - see Theorem 3b in [7j - entails then 

(2.16) P[T, > t] ~ cx'^TP 

for such processes, where c is an explicit constant. Notice that there are other Levy processes 
enjoying the property (12.151) . like subordinate stable processes, for which no refinement of 
(I2.14p seems available in the literature. In [4j, the precise estimate 

was obtained for every centered Levy processes with finite variance and every x > 0. 

We conclude this paragraph with spectrally negative Levy processes, where our survival 
analysis amounts to the study of a one-dimensional probability as in the discrete framework, 
and does not really require (I2.12p . Indeed, the passage-time process {T^jX > 0} is then 
a subordinator with infinite lifetime if Z does not drift to — oo, whose Laplace exponent 
$(A) = — logE[e~'*'"^^] is characterized by the law of Zi - see Chapter VII in [5j or Chapter 

9 in [20] for these basic facts. In particular, if $(A) ~ A''/(A) as A — )■ for some p G (0, 1) 
and /(A) some slowly varying function, then Theorem XIII. 5. 4. in [22] yields 

p[r. > t] ~ Y^^iit-')r'- 

Actually, the above condition on $ is equivalent to (I2.13P - see Proposition VII. 6 in [5J 
and notice that then necessarily p > 1/2. Hence, the above estimate is just a consequence 
of Rogozin's theorem with an explicit constant c^. Spectrally negative Levy processes also 
enjoy the following peculiar property, which follows easily from the Baxter-Donsker formula 
expressed with the exponent $: 



(2.17) Mr^ ^ Zr^ I Zr^ > 0. 

In particular, one has Mt = Zt\ Zt > for every t > ( 

depend on t (which is true only for the (l/p)-stable process). This latter identity which can 



In particular, one has Mt = Zt\ Zt > for every t > if f[Xt > 0] = p > 1/2 does not 
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be shown in many different ways - see Section 8 in [8] and the references therein, recovers 
the estimate (12.161) in this particular case: one has 

for some exphcit constant c > which is a/2 for the standard Brownian motion. We finally 
stress that (12.171) can be useful for spectrally negative processes such that f[Xt > 0] — )■ 1. 
For example, if X is an a— stable process with positive drift then (12.171) shows after a simple 
analysis that 

for some explicit c > 0, an estimate which cannot be obtained directly neither from Rogozin's 
theorem nor from the behavior of $ at zero, and which is also coherent with (12. 7p since 
P[Xt < 0] ~ ct^-". 

3. Recent advances 

3.1. Integrated random walks. An integrated random walk is the sequence of partial sums 

An = Si + - ■ ■+Sn, where {S'n, n > 1} is a random walk. As above we write Sn = Xi + - ■ ■+X„ 
and we denote by /i the law of the increment Xi = Si. Let 

Tx = inf {n > 1, An > x} 

be the first-passage time above x > 0. Since {Si < 0, . . . , 5*^ < 0} C {Ai < 0, . . . , A„ < 0} , 
the discussion made in Paragraph 12.11 show that P[T^ = +oo] > P[To = +oo] > as soon 
as Sn drifts to — oo. If fi is concentrated on M+, then An has non-negative increments and 

since An = nXi + {n — 1)X2- ■ ■ + X„ > {n/2)S[n/2]^ one has 

P[T, >n] = ¥[An <x] < ^nSin/2] < 2x] < e'^+niogne-^"^'^''^) ^ 

which shows that ¥\Tx > n] tends to zero superexponentially fast, unless /i is degenerate. If 
/i is not concentrated on M+, then choosing e such that /i(— oo, —e) > and using the same 
argument as in Paragraph 12. II entail P[Tx > ''^] > Pi^o > f^] ^ (Ai(— oo, — £:))'^P[Ta, > n] for 
every x > as soon as ke > x, so that 

P[r, > n] X P[ro > n] 

for every x > 0. 

Henceforth we will suppose that Sn does not drift to — oo and that fi is not concentrated 
on M"^. This entails that P[S'„ > 0] > and F[Sn < 0] > for every n > 1. We are interested 
in the rate of decay of P[To > n] to zero, and we will see that far much less is known than 
for random walks. 

The case of integrated simple random walks was first considered by Sinai [52], who showed 
the following 

Theorem (Sinai). Suppose that fi{+l} = fi{ — l} = 1/2. Then 

(3.1) P[To>n] X n-^/\ 
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The main idea lying behind this theorem is very simple. Let {rfc, A; > 0} be the a.s. finite 
sequence of return times of Sn to zero viz. tq = and r^ = mi{n > t^-i, Sn = 0}. On the 
one hand, the simplicity assumption entails that {Ar„,n > 1} is an integer- valued symmetric 
random walk, and that the a.s. identification 

{Ar,<0,...,Ar„<0} = {A<0,\fi = l... r4 

holds. The series X]n>i ^^^^[^-rn = 0] is convergent, and Sparre- Andersen's formula com- 
bined with the Tauberian theorem for monotonic sequences shows that 

> ^n] ~ 

■n 



On the other hand, the sequence r„ grows like n^ at infinity (more precisely, n~'^Tn converges 
in law to some positive (1/2)— stable law), so that after some residual analysis on the bivari- 
ate random walk {Ar^, r„} , one obtains the desired result. 

Since then several authors have tried to extend the validity of the estimate (13.11) to more 
general random walks, and the following conjecture is made in [H] and 



Conjecture 1. Suppose that /j, has finite second moment and zero mean. Then 

P[To > n] X n-i/^ 

A general result in this direction was obtained recently by Dembo and Gao ^16j, who 
proved the following 

Theorem (Dembo-Gao). Suppose that fi has finite second moment and zero mean. Then 
there exists an explicit constant K such that 

P[To >n] < Kn-^'^, n>l. 

If in addition 3/3 > such that /i(— oo, — r) > e^^^ for some r > and /x(— oo, —t) < e~^* 
as t ^ +00, then 

P[To > n] X n-^/\ 



The method used [16] is completely different from Sinai's and relies on a decomposition 
of the integrated walk at its supremum, which is somehow reminiscent of Sparre Andersen's 
argument, and too involved to be discussed here in detail. This method also allows for an 
elementary proof of the formula (12.41) in the symmetric and absolutely continuous case - see 
Proposition 1.4 therein. Notice that the result in [16] is stated in a different manner than 
ours - see Theorem 1.2 and Remark 1.3 therein - and that the assumption for the lower 
bound is slightly less stringent. But as in our formulation, this assumption essentially means 
that E[e~''"^i] < +oo for some /3 > and that the left tail fi{—oo, —x) has some regularity in 
the neighbourhood of zero - see Remark 1.1 in [16]. Shortly before this result, it was shown 
in [3] in using strong approximation that 

n-i/^(logn)-^ < F[To>n] < ^-^/^(logn)^ 

at infinity when Efe^l^^^l] < +oo for some /3 > 0. 
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Sinai's method of decomposing the integrated path along the excursions away from zero 
of the underlying random walk can be extended under some extra assumptions on the law of 
—Xi conditioned to be positive, which we denote by /i~. This idea had been used in [55] to 
show the same result as [16] in several situations like double-sided exponential, double-sided 
geometric, left- continuous, or lazy simple random walks, which all satisfy Dembo and Gao's 
assumption for the lower bound. Recently in [56] , this technique is combined with local limit 
theorems to obtain the following more precise result. 

Theorem (Vysotsky). Suppose that fi_ is exponential or that fi is left- continuous. If fi has 
finite variance and zero mean then there exists c = c{fi) > such that 

P[To > n] ~ cn-^/^. 

Up to now. Conjecture 1 is still unsolved in general, and the main challenge is to get rid of 
the extra assumptions on /i_. As can be seen from the central limit theorem, random walks 
with zero mean and finite variance are such that W[Sn > 0] — )■ 1/2. In view of the discussion 
made in Paragraph 12.11 it is hence natural to raise the more general 

Conjecture 2. Suppose that P[S'„ > 0] -)■ 1/2. Then 

¥[To>n] = n-^/^+°^^\ 

In this formulation, no assumption is made on the moments of /i and this enhances sharply 
the difficulty of the problem, which probably requires more combinatorial tools than the one 
used in the above references. 

We now turn to some situations where the persistence exponent of integrated random 
walks is not 1/4. We denote by T>{a) the set of probability measures attracted to some 
strictly a— stable law with a ^ 2 and we refer e.g. to Chapter XVII. 5 in [22] for more on the 
subject. Recall that if /i G T>{a), then E[|Xi|''] < oo for every s G [0, a) and that E[Xi] = 
if a > 1. The following is a consequence of the main result in [T6] : 



Theorem (Dembo-Gao). Suppose that fi G T>{a) for some a G (1,2). Then there exists an 
explicit constant K such that 

P[To >n] < /s:r^-(^~^/")/^ n>l. 

If in addition fi is attracted to a spectrally positive a— stable law, then 

P[To > n] X ^-(i-i/")/2_ 

The result in [16] is formulated in a different manner and is actually more general. In the 
case when fi G T>{a) however, the additional assumption made therein for the lower bound 
is equivalent to the spectral positivity of the attracting law - see e.g. Theorem XVII. 5.1 in 
^22j . With the help of Sinai's method, the estimate was recently refined in [5^ : 

Theorem (Vysotsky). Suppose that fi^ is exponential or that fi is left- continuous. If fi is 
normally attracted to some spectrally positive a— stable law with a G (1,2) then there exists 
c = c{fj,) > such that 

ro>n] ~ cn-(i-i/-)/2. 



12 F. AURZADA AND T. SIMON 

The result in [56] is stated differently but we see from [10] p. 382 that it is actually 
equivalent to our formulation. In view of our final discussion made in Paragraph 12.1^ it is 
very surprising that the case where /i is right-continuous and normally attracted to some 
spectrally negative a— stable law seems to be more difficult to handle than the dual situation 
where fi is left- continuous and normally attracted to some spectrally positive a— stable law. 

It is natural to ask what the persistence exponent should be when the limit stable law has 
negative jumps. Let us hence denote by V^a, p) the set of probabihty measures attracted to 
a strictly a— stable law with positivity parameter p. Notice - see e.g. [58] for details and recall 
that we excluded the one-sided case - that p G (0, 1) for a G (0, 1], that p G [1 — 1/a, 1/a] 
for a G (1,2), and that spectrally positive a— stable laws with a G (1,2) are such that 
p = 1 — 1/a. Besides, one has P[S'„ > 0] — > p whenever p G V{a,p) - see e.g. Theorem 
XVII. 5.1 in j22]. In view of the discussion made in Paragraph 12. H it is natural to raise the 
general 

Conjecture 3. Suppose that P[S'„ > 0] -)■ p G (0, 1). Then 

¥[To >n] = n-^/2+o(i)_ 

A weak version of the conjecture is to get the asymptotic under the assumption that 
p G V{a,p), just like Conjecture 1 which deals with the case a = 2 and p = 1/2. Under 
the stronger assumption that p is normally attracted, one may also wonder if a more precise 
behavior could not be obtained, as in 



3.2. Integrated Levy processes. In this section we consider the process 

At = f Z, ds, t > 0, 
Jo 

where {Zt, t > 0} is a real Levy process starting from zero. We set 

T, = inf {t >0, At> x} = inf {t > 0, At = x} 

for its first passage time above x > 0. Contrary to f l2.10p . there is no straightforward bound 
between F[Tx > t] and an analogous quantity involving some iterated random walk, so that 
the results of the previous paragraph cannot be used directly. 

The process [Z, A) is strongly Markovian and we set P(2,a) for its law starting from {z, a), 
with the simplified notation P = P(o,o)- It is clear by the right-continuity of Z that 

P{2,o)[^o = 0] = or 1 according as z > or 2; < 0. 

One has also P[To = 0] = or 1 by the 0-1 law, but to obtain a criterion for the regularity 
of the upper half-plane for {Z, A) is an open problem which does not seem obvious. Since 
integrated Levy process all have finite variation, one might think that this criterion will be 
different from the aforementioned Theorem 22 in [20j. 

If Z drifts to +00, then it is clear that T^ < +00 a.s. for every x > 0. On the other 
hand, if Z drifts to +00 then its last passage time above zero can be made arbitrarily small 
so that one will have P[Tj, = +00] > for every a; > 0. When Z oscillates then probably 
one has P[Ti. = +00] = for every x > 0, but there is not a direct answer to this question. 
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In general there is no result of basic fluctuation theory available for integrated Levy processes. 

In this paragraph we will consider two examples where the persistence exponent can be 
computed. The first one is the integrated Brownian motion and originates from Kolmogoroff 
[35] , in relation with the two-dimensional generator 

1 92 d 

2 dx"^ dy 

and the associated Fokker-Planck equation. Notice that [35j actually deals with the more 

general n-times integrated Brownian motion. The process {B, A) is a Gaussian Markov 

process whose transition density can be computed via the covariance matrix: under P = P(o,o) 

one has 

\/3 
pt{b,a) = -^exp[-2b'^/t + 6ab/t'^ -6a^/t% 

Tit 

and the expression under P(x,2/) follows by translation. By the (3/2)-self-similarity of A, one 
has 

under P, so that our persistence problem amounts to find the asymptotic of P[Ti > t], a 
question which dates back to Uhlenbeck and Wang in 1945. Notice that the above identity 
also yields F[Tq = 0] = 1. Among other formula, the following was obtained by McKean in 
an analytical way - see (3.1) in [41]: 

F(-i,o) [To G dt^Br, G dx\ = _^^^-i2/m-.+.^) / e-'^/'y-'/'dy l^.>oydtdx. 

nV^Tct^ \ Jn / 



This formula is the key argument to the following result which is proved separately in [23] 
and [2H]: 

Theorem (Goldman, Isozaki- Watanabe) . With the above notations there exists c > such 
that 

(3.2) P[Ti >t] ~ cr^/^ 

Notice that by self-similarity this result has the more general formulation 

F[As < X, Vs G [0,t]] = P[T, > t] ~ cx^/^t"^/^ as t/x^/^ -> +00. 

In particular, one has 

F[As < X, Vs G [0, 1]] ~ cx^/^ as X ^ 0, 

which is a lower tail probability statement as mentioned in the introduction. By strong 
approximation and as a consequence of his result on integrated simple random walks, Sinai 
had obtained in [52J the rougher estimate 

p[ri > t] X t-^/^. 
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The argument in [23] is analytical and relies on integral equations. It also yields a complicated 
explicit expression for the law of Ti : 

P [T, E dt] = P(o,-i) [To G dt] = (^4^,^-'^"'^ 



9 AT pco pco pt \ 

+ ^— / / / P(-.,o) [t-ToE ds, Bto e dx] e-6/«'-2-'/«' smh{6z/ s^)zdzds dt, 

T^ Jo Jo Jo ' J 

where the quantity P(_2.o) [t — Tq E ds, B^q E dx] can be expressed via McKean's above for- 
mula and self-similarity. Notice that both McKean and Goldman's formulae have been 
generalized by Lachal [3S], who obtained an explicit formula for P(fe,a) [To E dt, Bt^ E dx] for 
any {b,a). The asymptotic analysis which is carried out in Proposition 2 of |23] gives the 
right speed of convergence for P[T E dt], but the resulting value of the constant c in (13. 2 p is 
erroneous because of some inaccuracies in the change of variable. The right value is 

, _ 3^/^r(2/3) 

7r2i3/i2r(3/4) ^ 

and follows also from the simpler probabilistic method of |28], which relies on the Markov 
property and a Tauberian argument. This method yields the more general estimate 

with some explicit Cb^^ > - see (1.12) in |2B]- It also allows to handle first-passage time 
asymptotics for fluctuating homogeneous additive functionals of Brownian motion, with a 
persistence exponent depending smoothly on the skewness of the functional - see Corollary 
1 in [27]. We finally notice that the estimate (13. 2 p was also obtained (with an erroneous 
constant c) in [13] after solving some Krein-Kramers differential equation in the context of 
semifiexible polymers in the half-plane. This latter method was generalised in pJJJ to give 
another computation of the persistence exponent for fluctuating homogeneous additive func- 
tionals of Brownian motion, in the context of survival of a diffusing particle in a transverse 
shear flow. 

We now turn to integrated strictly a-stable Levy processes, which form the natural gen- 
eralisation of integrated Brownian motion. The bivariate process {Z, A) is then a stable 
Markov process, where the stability property has the general meaning which is given in the 
monograph [l6]. In particular, one can check from the Levy-Khintchine formula - see e.g. 
Proposition 3.4.1 in [16] - that 

A, = {l + a)-^Zi 

for every t > 0, so that Ai is a a— stable variable with the same positivity parameter as 
Zi. On the other hand there does not seem to exist an explicit formula for the density of 
the bivariate random variable {Zi,Ai) except in the Gaussian case a = 2. The univariate 
process A is (1 + 1/a)— self-similar, so that with the above notations one has 
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under P and we need only to study the asymptotic behaviour of P[Ti > t]. Notice that again 
this identity yields P[To = 0] = 1. As a consequence of the main result of |50j, the persistence 
exponent of A can be computed within a specific sub-class: 

Theorem (Simon). Let {Zt, t > 0} be a strictly a-stahle Levy process with a G (1, 2). With 
the above notations, there exists a positive constant K such that 

p[Ti >t] < fs:t-a-V")/2^ ^>o. 

// in addition Z is spectrally positive, then 

P[Ti>t] = t-(l-l/")/2+o{l)_ 

The main result in [SO] deals with more general homogeneous functionals of stable Levy 
processes, extending the results of [26]. It also provides some explicit lower bound with a 
logarithmic term, which entails the following criterion for the finiteness of fractional moments 
of Ti in the spectrally positive case: 

W(ri\ < cx) ^ -(« + 1) < s < (a - l)/2a. 

The method of [50] is an adaptation of Sinai's in continuous time, relying on the bivariate 
Levy process {(rt,y4T-J, t > 0} with {rt, t > 0} the inverse local time at zero (which exists 
because a > 1). The upper bound relies on the Wiener-Hopf factorization method as in [26], 
and the fact that A^ is a symmetric (a — l)/(a + l)-stable variable whatever Ai's positivity 
parameter is. For the lower bound, the crucial fact is that a.s. 

{Ar, < 1, Vs G [0,t]} = {A, < 1, Vs G [0,rj}, 

which allows to study the probability of the right event with the help of (I2.16p . This latter 
identity is true only in the spectrally positive case. 

It is a natural question to find the persistence exponent for all integrated stable Levy 
processes. If Z is an a-stable subordinator, then A is an increasing process and one has 

where the estimate follows e.g. from (2.4.30) in [25]. This superexponential speed is coherent 
with the previous discussion for integrated one-sided random walks. If —Z is an a-stable 
subordinator, then T^ = +oo a.s. for every a; > 0. The general problem remains open, and 
in view of Bingham's estimate (I2.16P it is natural to raise the precise 

Conjecture 4. Let Z be a strictly a— stable Levy process such that P[2'i > 0] = p G (0, 1). 
Then there exists c > such that 



We see that the main result of [50] shows the conjecture in the spectrally positive case viz. 
a > 1 and p = 1 — 1/a, except the existence of the constant for which the refined methods 
of [SSj could perhaps be adapted. In the general case, a simple scenario with an initial 
downward intention is described in [TS] p. 3 to show that the speed must be larger than t~P. 
Unfortunately, the attempts made therein to obtain a better bound were unjustified and in 
[51] a method relying on decorrelation inequalities, inspired by jT^, is carried out to show 
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that the upper bound in the main result of [50] is not the optimal one when Z has negative 
jumps and self-similarity index close to 1: 

Theorem (Simon). For every ao, c > there exists k > such that for every a G [«0;2] 
and every strictly a— stable process Z with positivity parameter p E [c A (1 — l/a), l/a V 1] 
one has 

liminft'^PfTi > t] =0. 

At present, it does not seem that the different arguments developed in [181 EDI EI] can 
provide any clue to get the conjectured persistence exponent p/2. It is quite surprising that 
when Z is spectrally positive the persistence exponent is easier to get for A than for —A. 
Indeed, the simple identity (12.171) entails readily that the persistence exponent of —Z is 1/a, 
whereas it is not a straightforward task to show that the persistence exponent of Z is 1 — 1/ a. 

It seems plausible that a version without constant of Conjecture 4 could be proved in 
refining the results of [16] and using some discrete approximation. In view of Conjecture 3 
and (I2.14p . this leads us to the following general question on integrated Levy processes: 

Conjecture 5. Let Z he a Levy process such that F[Zt > 0] — ?■ p G (0, 1) as t -^ +cxd. Then 
for every x > 



3.3. Fractionally integrated Levy processes. In this section we consider processes of 
the type 

(3.3) Af = ^^^ ^ (i - sf dZ,, t > 0, 

where {Z^, t > 0} is a real Levy process starting from zero and /3 > — 1. The above convo- 
lution product makes sense for every (3 > 0, and can also be defined for some negative f3 
depending on the law of Z. If Z is strictly a— stable for instance, then A^ is well-defined 
for every (3 > —(1 A l/a) and is then a stable process in the sense of [46], with continuous 
paths iff a = 2 or /3 > and with a.s. locally unbounded paths iff a < 2 and /3 < - see 
Chapter 10 in [IB]- When Z is strictly a— stable, it is customary to write (i = H — 1/a with 
H > the so-called Hurst parameter, and A^ is then an if —self-similar process. One can 
view /3-fractionally integrated Levy processes as the natural generalization of n times inte- 
grated Levy processes, since an integration by parts shows that the latter form the subclass 
P = n. In this particular case the (n -|- l)-dimensional process [Z, A^, . . . , A") is Feller, but 
there is no multidimensional Markov property when /3 is not an integer because then the 
fractional integration takes the whole memory of the driving process into account. In the lit- 
erature fractionally integrated Levy processes are often called Riemann-Liouville processes, 
a denomination which is originally due to Levy. 

In the Brownian case Z = B, the process A'^ is closely connected to the fractional Brownian 
motion {B^ , t > O}, which we recall to be the centered Gaussian process with covariance 
function 

E [BfB^] = ^ (t^^ + s^^ -\t- s|2^) , t,s>0. 
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Fractional Brownian motion can be written as the independent sum 

(3.4) Sf = CH U"^-''' + r^it + 3)^-1/2 - s^-V2) d^^j , t > 0, 

with ch = {H2'^^Tx /V{H + l/2)r(l — H)y^'^ the normahzation constant, which shows that 
its paths are continuous a.s. The process B^ gives insight on the long-range increments of 
j^H-i/2 gjj^pg j^ pg^j^ i^g proved in analysing the covariance function that 



-CXD 



(3.5) {<-;/^ - A^-y^ t > 0} A {c-^'Bf, t > 0} , n ^ 

(with, of course, an equality in law for every u when H = 1/2). Recall that B^ is defined for 
H G (0, 1] only and that B^ is simply the linear function t i— )■ tN with A^ a standard normal 
variable. Fractional Brownian motion can be shown [46] to be the unique if —self-similar 
Gaussian process with stationary increments, whence its greater importance in modeling 
than Riemann-Liouville processes. Setting 

Tf = inf {t > 0, Bf >x} = inf {t > 0, 5f = x} 

for every x > 0, one has T^ = x^^^T^ by self-similarity, which shows that T^ = a.s. and 
that the survival analysis of B^ is reduced to the behavior of F[T^ > t] only. Among other 
results, the following was obtained in 



Theorem (Molchan). For every H G (0, 1] one has 

F[Tf >t] = t^-i+°(i). 

The main argument of [52], partly inspired by Brownian fluctuation theory, is to quantify 
the correlation between T/^, the last zero of B^ on [0, 1], the positive sojourn time of B^ 
on [0, 1], and the inverse exponential functional 

t^ J,^ = (^*e^?d.y\ 

whose asymptotic behavior in expectation can be precisely analysed. The general link be- 
tween J/^ and T^ is explained by the heuristical fact that if Tf^ > t then B^ has drifted 
towards — oo rather rapidly, so that J/^ is big. Conversely if T/^ < t then B^ has been 
close zero for a positive fraction of time, so that J/^ is small. However, the analysis of 
{E[J/^], t > 0} which is performed in [12] is very specific to the stationary increments of 
fractional Brownian motion, and does not seem to be suitable for Riemann-Liouville pro- 
cesses. 

It is clear that ¥[T^ > i] — t- 1/2 by the above remark on the case H = 1, and it is natural 
to raise the 

Conjecture 6. For every H G (0, 1) one has 

P[Tf > t] X t"-\ 
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In view of the previously stated precise results, one may also ask if P[T/^ > t] ~ ct^^^, 
but from empirical studies which were communicated to us by A. A. Novikov, it seems that 
such an exact behaviour does not hold in general. According to these studies, at least for 
H = 3/4 one should have 

liminf t^"^P[Tf > t] < limsupt^"^P[Tf > t]. 

The above conjecture was recently addressed in [Ij, where the following partial result is 
obtained by a refinement of Molchan's method: 

Theorem (Aurzada). There exists c > such that 

{\ogt)-H^-^ < P[Tf >t] < {\ogtyt^-\ t^+oo. 

We now turn to first passage asymptotics for a class of Gaussian stationary processes 
(GSP) which are related to the persistence of fractionally integrated Brownian motion. In 
[38j the Lamperti transformation 

(3.6) Lf = e-*^S5, teM, 

a centered GSP, is studied in connection with the persistence exponent of B^ and it is shown 
that 

logP [Lf < 0, Vs G [0,t]] = logP [5f < 0, Vs G [l,e*]] ~ t{H -1). 

Since Molchan's theorem means that logP [B^ < 1, Vs G [0,e*]] ~ t{H — 1), one simply 
needs to switch the and the 1 in the latter probability to obtain this result, which is 
justified in p8] through a refined use of Slepian's lemma. The Lamperti transformation is 
also a fruitful method in the reverse direction, and this was observed in [3J to investigate the 
persistence exponents of Riemann-Liouville processes. Set Z = B in (13. 3p and introduce the 
notation 

T^ = inf h > 0, Af > x} = inf It > 0, Af = x\ 

for every f3,x > 0. Notice that T^ = x^/^^"'"^^^Tf by self-similarity, which shows that Tq = 
a.s. and that the survival analysis of B^ is reduced to the behavior of P[Tf > t] only. The 
process 

(3.7) Y/" = e-^^Aj-'/', teR, 

a centered GSP for every H > which reduces to the stationary Ornstein-Uhlenbeck process 
when H = 1/2, plays a key role in the following 

Theorem (Aurzada-Dereich). There exists a non-increasing function /3 h-)> 6 {(3) such that 

for every (3 > 0. Besides one has 9{oo) = — lim^^oo t~^ logP [F^ < 0, Vs G [0, t]] , where Y is 
the centered GSP with correlation E[Yo^t] = l/cosh(t/2). 
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The important process Y is mentioned in [39] in the context of diffusion equation with 
white noise initial conditions, and in [T7] in relation with the positivity of random polynomials 
with large even degree. More details will be given in the next section. In p8], the process 
Y is also viewed as the Lamperti transformation of the curious smooth (l/2)-self-similar 
Gaussian process 



Xt = V2t^ / BuC-'^^du, t > 0, 
Jo 



which shares the same time inversion property {Xf, t > 0} = {t ^X^-i, t > 0} as Brownian 
motion. In [3] only the upper bound 

(3.8) - lim t-^ log P [y, < 0, V s G [0, t]] < 9{oo) 



t—^ca 



is proved via Slepian's lemma, but one can easily show that Y^ — )■ F as if — ?■ oo in analysing 
the covariance function, so that the inequality in (13. 8p is actually an equality, the limit on 
the left-hand side being a supremum. 

The above result has several interesting consequences. First, it shows that the persistence 
exponent is a non- increasing function of the order of integration for Brownian Riemann- 
Liouville processes. The fact that smoother processes have more probability to survive is 
believed to be a kind of universal feature. Second, it entails that 6'(/3) > 6{1) = 1/4 for 
every f3 < 1, so that by Molchan's result the persistence exponents of B^ and A^^^^'^ do 
not coincide whenever H > 3/4. We actually believe that 6{H — 1/2) > 1 — H ioi every 
H e (1/2, 1) and some reasons for that will be given soon afterwards. Last, it entails that 
6{oo) < 6{1) = 0.25, which improves the bound ^(oo) < 0.325 obtained in [3S] via another 
Slepian's inequality. In [39] the numerical value 9{oo) ~ 0.1875 is suggested, whereas in 
[T7] the value ^(oo) = 0.19 ± 0.01 is obtained by simulations. It is a tantalizing question to 
compute the function 9{/3) for every positive /3 ^ {0, 1}, as well as its limit 6^(00). The lower 
bound ^(00) > 0.125 is obtained in [37] with the help of a certain Gaussian comparison in- 
equality, and in [H] this lower bound is improved into 9{oo) > l/4v3 > 0.144, in comparing 
Y with a linear time-change of the so-called Wong process [57J, which is the GSP associated 
with integrated Brownian motion. 

Let us give some more details on the correlation function Cnit) = K[Yf^Y/^] of the Lam- 
perti transform introduced in (13. 7p . It is given by 



Jo 



for every t > 0, and a simple analysis shows that for every H G (0, 1] and t > one has 



Cnit) < cosh{Ht) - 2^^-i(sinh(t/2))^^ = E[L^ L 



t J' 



where L^ is the Lamperti transform defined in (13. 6p . By Slepian's lemma and Molchan's 
theorem, this entails 

(3.9) e{H-l/2) > 1-H 
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for every H G [1/2, 1], with an equality ii H = 1/2. From the above we know that (13.91) is 
strict if if > 3/4. It is very hkely that 6{H - 1/2) exists for every H G (0, 1/2) and that 
( 13. 9p is strict for every H ^ 1/2. In view of the above discussion, it is natural to raise the 

Conjecture 7. The function /3 h-)> 9{I3) is convex decreasing. 

It is also interesting to look at the first-order expansion at zero 



2H 



Til-H)Til/2 + H) ft\ , ^,^2H^ 



which should be compared with 

2H 



E[L^Lf] = 1 - (0 + o{t^^] 



One might actually wonder if, for every H G (0, 1/2) U (1/2, 1), the sole inequality 

Til-H)Til/2 + H) 

V^ 
could not be enough to prove that 9{H — 1/2) > 1 — H. Indeed, even though the asymptotic 
analysis of large excursions of continuous GSP's is known to be a hard problem in general, 
which should involve the whole correlation structure, there are also reasons to believe that 
its behavior at zero plays a prominent role when comparing self-similar processes having the 
same sample path regularity. For H = 1 one has 

Ci(t) = 1 + ^logt + ait') 

and from this logarithmic behaviour, which is related to the local growth of Brownian motion, 
one might view half-integrated Brownian motion as a kind of boundary object, analogous to 
the Cauchy process among stable Levy processes, and where the survival analysis should be 
different than in the other cases. This boundary phenomenon has actually been observed by 
physicists for a while, and we will give more details in the next section. For H > 1 one has 

which entails by the Ito-Rice formula [29] that the zero-crossings of the process Y^ are 
isolated points, whose number N^ on an interval of length t has expectation 



E[<] ' ' ^ 



271 \j {H- 1) 

The interesting fact that H i— t- E[A^/^] is decreasing might well be related to the non- 
increasingness oi H \-^ 9{H — 1/2) = — limi^_|_oo ^~^ logP[A^/^ = 0]. Although no mathe- 
matical result seems available, all empirical results in the physical literature [391 HH] show 
namely that for smooth centered GSP's with positive correlation function F and zero-crossing 
number Nt, the greater — F"(0), the greater 9 = — limi^_|_oo t~^ logP[A'"i = 0]. 
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Let us shortly notice that integrated Brownian motion is somehow critical as far as the 
small values of the length of the excursions of the associated Lamperti process are concerned. 
As a consequence of the Longuet-Higgins formula [57] and a second-order expansion of Cnit) 
at zero, one has indeed 



-■5 



+CX) if if < 3/2 
limr^P[iVf ^0] = <J 37/96 ii H = 3/2 

if if > 3/2. 



t^o 



However, this has probably less to deal with the fact that 6{H — 1/2) is computable only for 
H = 3/2, since analysing small and large values of the length of the excursions of continuous 
GSP's are very different problems. 

We now go back to general fractionally integrated Levy processes and state the following 
result of ^, which is obtained by strong approximation and drift transformations of Gaussian 
processes: 

Theorem ( Aurzada-Dereich) . For every P > 0, the persistence exponent 6{P) is the same 
among all (3— fractionally integrated Levy processes such that t t-)- E[e*'^^] is defined in an 
open neighbourhood of zero. 

The main result of [3], to which we refer for details, is more general and handles frac- 
tionally integrated random walks as well as more general Volterra kernels. The finiteness of 
exponential moments play an important role therein due to the use of strong approximation, 
but in view of the aforementioned results and conjectures on integrated random walks and 
Levy processes, it is natural to raise the 

Conjecture 8. For every /3 > 0, the exponent 6{P) is the same among all P— fractionally 
integrated Levy processes with finite variance. 

We conclude this paragraph with the integrated fractional Brownian motion 

If = [ Bfds, t>0. 
Jo 

This is a centered Gaussian process with positive correlation function, and with the help of 
its Lamperti transformation it is easily shown that for every H G (0, 1) there exists p{H) > 
such that 

p[/f < 1, VsG [o,t]] = r^^^^+^w. 

Numerical simulations [15] suggest the following: 



Conjecture 9 (Khokhlov-Molchan). One has p{H) = if (1 — H). 

This expected value, which is symmetric with respect to if = 1/2, is very surprising 
because it is known that fractional Brownian motions with Hurst index smaller or greater 
than 1/2 are very different processes from several viewpoints. This does not match either 
the above heuristic discussion on smooth GSP's since the correlation function F^ of the 
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Lamperti process of I^ has a first-order expansion at zero 

F^it) = 1 - (±lpl + o{t^). 

In particular, the number N^^ of zero-crossings of this Lamperti process has an expectation 
E[A^/^] = ty/l — H'^/n which decreases with H, and one might think that p{H) also decreases. 
However integrated fractional Brownian motion may well be more the exception than the rule 
for this kind of questions, because of its complicated correlation structure. In |15] it is argued 
that the difference between H < 1/2 and H > 1/2 should be observed at the logarithmic 
level and it is shown in pS], with a detailed analysis, that cH{l — H) < p{H) < 1 — H for 
some c G (0, 1). Other bounds such as p{H) > p{l — H) for every H < 1/2 were also recently 
presented in [H] (with numerical explanations), neither proving nor disproving the above 
conjecture. 

3.4. Other processes. Integrated random walks which were considered previously can be 
written as the weighted sum 

n 

^n = ^(n+l-2)Xi, 

i=l 

where {Xj, z > 1} is an i.i.d. sequence. The above weights depend on i and n so that the 
increments of An are neither stationary nor independent. It is natural to consider persistence 
problem for other weighted sums like 



y^o'»^i 



where {cjj, z > 1} is a deterministic sequence and {Xj, i > 1} an i.i.d. family. The situation 
is a bit simpler than for integrated random walks because the non-stationary increments of 
E are independent, nevertheless the sequence {ai, i > 1} has also more generality. Setting 
p for the law of Xi and introducing 

T, = inf{n>l, S„>0}, 

the following has been obtained in |2] via strong approximation techniques: 

Theorem (Aurzada-Baumgarten). Suppose that p is centered and that its Laplace transform 
is finite in an open neighbourhood of zero. Suppose that {a^, i > 1} is increasing with an ^ n^ 
for some p > 0. Then 

F[T^>n] = n-(P+i/2)+o{i)_ 

The results of [2] are more precise and allow other weight functions not necessarily in- 
creasing when p is Gaussian. The general case reduces to the Gaussian one via strong 
approximation and drift transformations, as in |3]. A universal speed not depending on p 
can also be obtained for weight functions growing faster than polynomials, like e"^ for some 
7 < 1/4, and the surviving probability is then exponentially small. However, it is also shown 
in [2] that the speed does depend on p for weight functions growing too fast, like e"'"' for 
some 7 > 1. It would be interesting to find a's critical growth rate for the universality of the 
speed. The following is also a natural question. 
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Conjecture 10. Suppose that fi is centered and has finite variance. Suppose that an ^ n^ 
for some p > 0. Then 

Let us conclude this paragraph with iterated Levy processes, which are processes of the 
type {X o |l^t|, t > 0} with X, Y two independent real Levy processes starting from zero. If 
|y| is a subordinator, then X o |y| is another Levy process which is called a subordinate 
Levy process, a notion introduced by Bochner in the context of harmonic analysis. Iterated 
Levy processes were introduced by Burdzy in the Brownian framework and can be viewed as 
a generalisation of subordinate Levy processes. They are known to have strong connections 
with PDE's of higher order, especially through their first passage times. Let us introduce 

T, = inf {t > 0, Xo ir^l > 1} 

for all a; > 0, and notice that as for integrated Levy processes the law of T^. is difficult to 
study in general since X o |y| is non-Markov. Among other results, the following has been 
obtained in |1]: 

Theorem (Baumgarten). Suppose that the random variables |Xi|" and \Yi\" have exponen- 
tial moments for some a > and that E[Xi] = 0. Then 

with e = 1/4 ifE[Yi] = ande = 1/2 ifE[Yi] ^ 0. 

In particular, the persistence exponent of iterated Brownian motion is 1/4. The strong 
dichotomy between the situations where Yi is centered and non-centered is not surprising 
since in the former case \Yt\ grows roughly like y/t whereas in the latter case \Yt\ grows like t. 
As in many above statements, the question of replacing the exponential moment condition 
by the sole assumption of finite variance remains open. 

4. Some connections with physics 

4.1. Regular points of inviscid Burgers equation with self-similar initial data. The 

statistical study of the one-dimensional Burgers equation 

(4.1) dtu + udxU = vdxxU 

with viscosity z/ > and an initial condition uo{x) := m(0, x) = X^ given by a self-similar 
stochastic process {X^, x G M} has been initiated in the papers |l9l |53]. Though this equa- 
tion is accorded to be an unrealistic physical model for turbulence, the competition between 
the irregularities of X and the irregularities generated by (14. ip remains an interesting math- 
ematical study. In the inviscid limit z/ = 0, the Hopf-Cole solution to (14.11) is given by 

u(t,x) = '^-^'^^'^'^ 

for every t > 0,z G M, where a(t, x) = max{|/ G M, Ct{y) < xt~^} and C is the right- 
derivative of the convex minorant of the function 



ry 
Ft-.y ^ / (X, + xt-^) dx. 
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This variational formula is obtained in considering the explicit solution to (14. ip which can 
be obtained ioi u > and letting z/ — )■ - see [IHl El EB IS3] for details. Notice that a{t, x) 
is well-defined only if 

(4.2) l^l'^Xj; — )■ a.s. when x — )■ ±oo. 

The function x ^ a{t, x) is right-continuous but not continuous in general, and the so-called 
Lagrangian regular points at time t > are defined as the set 



Ct = {a(t,a;), x G M and a{t,x—) = a{t^x)} 

which consists of the points where Ft equals its convex minorant. In physical terms, the set 
Ct describes the initial locations of the particles which have not been shocked up to time 
t. It is easy to see that when X is a self-similar process, the map t ^ Ct has also some 
self-similarity which makes the a.s. Hausdorff dimension of Ct independent of t > 0. Setting 
C = C\ and "Dim" for "Hausdorff dimension", the following is stated in \0\: 

Conjecture 11 (Aurell-Frisch-She). Suppose that X is the fractional Brownian motion with 
Hurst parameter // G (0, 1). Then Dim C = H a.s. 

Notice that in the above, the fractional Brownian motion is defined over the whole M, and 
coincides with the two-sided Brownian motion when H = 1/2. This conjecture remains open 
in general, but the following has been shown: 

Theorem (Handa-Sinai, Bertoin). Suppose that X is the two-sided Brownian motion. Then 
Dim C =1/2 a.s. 

This result was first stated in [53] , although no strict proof is given therein for the lower 
bound Dim C > 1/2 a.s. A simple argument based on integration by parts and Frostman's 
lemma is presented in [2l], which yields the general lower bound Dim C > if a.s. when X 
is the fractional Brownian motion with Hurst parameter H. In [6j, the exact computation 
of the Hausdorff dimension follows as a simple corollary to the more general result that 
X I—)- a{l,x) has stationary and independent increments with explicit Laplace transform. 
This result extends to Levy processes with no positive jumps satisfying the growth condition 
(14. 2p . In particular one has the 

Theorem (Bertoin). Suppose that X is a two-sided a -stable spectrally negative Levy process 
with index a G (1,2). Then Dim C = 1/a a.s. 

We now briefly describe the link between an upper bound for Dim C and the computation 
of certain persistence exponents, in the self-similar framework. This is the original argument 
of [53] for Brownian motion and it extends to fractional Brownian motion [l5] or a-stable 
Levy processes [51j- Speciflcally, it can be shown by the Borel-Cantelli lemma and some 
elementary inequalities that if 

r(X, + x)dx > -(5^+^, Vy G [-1, 1] < 5i-^+°(i) 

'0 

as (5 —7- 0, where H is X's self-similarity index and K G [0, 1], then Dim C < K a.s. If 
X = B^ is the fractional Brownian motion, then the drift appearing in (14. 3 p can be removed 



(4.3) P 
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by quasi-invariance and some analysis, and one sees by symmetry and self-similarity that 
the required estimate to get the upper bound in Conjecture 11 is 



P 



r B^dx < 1, yye[-t,t] 
Jo 



< t'-H+oH) ^^+oO. 



This is a "two-sided" persistence problem, and the above estimate is formulated as a conjec- 
ture in |l5l Us] , with an equality instead of the inequality. Notice that this latter problem is 
independent of Conjecture 9 since the increments of B^ are correlated. Actually, even the 
sole existence of the persistence exponent for integrated double-sided fractional Brownian 
motion has not yet been established. 

If X is a two-sided a-stable Levy process with a G (1,2), some tedious analysis shows that 
the drift appearing in (14. 3 p can also be removed (this would not be the case for a G (0, 1]). By 
self- similarity and independence of the positive and negative increments of Z the inequality 
(14. 3 p amounts then to 

P[Ti > t] < t(i-^)/2+"(i) 

at infinity, where Ti is the first-passage time at 1 of the integral oi Z = —Z. In particular the 
validity of Conjecture 4 would lead to Dim C < p a.s. where p = P[Zi > 0] is the positivity 
parameter of Z. In view of the above theorem, it is natural to state the 

Conjecture 12. Suppose that X is a two-sided a-stable Levy process with a G (1,2) and 
positivity parameter p G [1 — l/a, 1/a]. Then Dim C = p a.s. 

This Hausdorff dimension depending on the positivity parameter and not on the self- 
similarity index is different from the value 1/a which had been conjectured in [31] through 
multifractal analysis. The invalidity of the latter when a is close to 1 was proved in [51] with 
the help of the above Sinai's approach: 

Theorem (Simon). For every c < 1 there exists ao > 1 such that for every a G (l,ao) 
and every p G [1 — l/a,c A 1/a], if X is a two-sided a-stable Levy process with positivity 
parameter p, then Dim C < 1/a a.s. 

We conclude this paragraph in mentioning that the lower bound in Conjecture 12, which 
should be obtained by different and yet unknown arguments, would lead by comparison to 
the lower bound of Conjecture 4 in the case a > 1. 

4.2. Positivity of random polynomials and diffusion equation. A classical question 
dating back to the beginning of probability theory is to understand the distribution of the 
roots of random polynomials. Consider 

n-l 

P„(X) = Y^^.x' 

1=0 

with large even degree where {^j, i > 0} is some i.i.d. sequence and X the deterministic 
variable, and set N^ for the number of its real roots. Among other results, the following was 
recently obtained in [T7] : 
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Theorem (Dembo-Poonen-Shao-Zeitouni). Suppose that ^i is centered and has polynomial 
moments of all order. Then 

(4.4) P[iV2„+i = 0] = n-^^+°(i) 

where b = — liiiii^oo^ ""^ logP [y^ < 0, Vs G [0,t]] , with Y the centered GSP with correlation 

E[YoYt] = l/cosh{t/2). 

In the above, the exact value of b is unknown and numerical simulations suggest 46 = 
0.79 ± 0.03 - see [IT]. It is remarkable that this constant b relates to n times integrated 
Brownian motion. We saw indeed in Paragraph 13.31 that b = lim„_^_|_oo^(^) where 6{n) is the 
persistence exponent of the process 

1 /■* 
t ^ — (t- s)"" dBs, t > 0. 

The problem of computing 6{n) for n > 1 is believed to be very challenging. Numerical 
simulations |40j suggest 9(2) = 0.231 ± 0.01. 

Let us give some insight on the proof of the above result. The hard part is to show (14. 4p 
when ^1 ~ Af{0, 1). The general case follows by strong approximation, whence the assumption 
made on the moments, but notice that it is also conjectured in [17] that (14. 4p should hold 
under the sole condition that ^i is centered and has finite variance. When ^i ~ A/'(0,1), 
the process x — )■ Pn{x), which is the so-called Kac's polynomial, is centered Gaussian with 
covariance E[P„(x)P„(|/)] = 1 + ■ ■ ■ + {xy)"'~^. Its correlation function is given by 



{xyr - 1 



xy — 1 



[x'-l){y'-l] 



(a;2"-l)(y2"-l) 



for every x,y ^ ±1. This function is invariant under the transformations {x,y) H- (— x, —y) 
and {x,y) h-)- (1/x, l/y), and an involved argument based on Slepian's lemma shows that 

p[iV2„+i = o] = (P[iv£iii = o])^+°« 

where Nn ' is the number of roots of P„ on [0, 1]. The link between A^2n'+i ^°^ large n and 
the zero-crossings of Y is established after changing the variable x = e~* and isolating the 
contributions for small t. The latter is the crucial step, since it follows from the singularities 
of the correlation function that the density of the real roots of P„ around ±1 is very big for 
large n. Notice that the link between Y and Nn is rather easily understood at the expectation 
level: a classical formula due to Kac - see PU] for details - yields 

E[iv£i] ~ ^logn, 

whereas the Ito- Rice's formula shows that if Nt is the number of zero-crossings of Y on [0, t], 
then 'E[Nt] ~ t/2Tr. The problem of evaluating P[A^2n+;^ = 0] is however much more intricate 

than the sole estimation of E[A^2nVM5 i^ analogy with what happens for the zero-crossings of 
Gaussian stationary processes. 
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It is also a remarkable fact that the above constant b appears as the persistence exponent 
of another, apparently disconnected random evolution phenomenon, which is studied in 
[391 ill ST]. Consider the heat equation on W^ 

(4,5) I^^A.. 

with random initial condition Ud{x, 0) = W{x) a d- dimensional white noise. Integrating 
along the heat kernel, it is easy to see that for every x G M'^ the solution t ^ Ud{x,t) to fl4.5p 
is a (— (i/4)-self-similar centered Gaussian process with covariance function 

E[udit,x)udis,x)] = ^^^^^^^^,/2 - 

In particular the law of {ud{t,x), t > 0} does not depend on x, which is also clear from 
the white noise initial condition. The Lamperti transformation t i-> (27r)^/^e*/^M(i(x, e*) is a 
centered GSP with correlation function l/(cosh(t/2))'^/^, and this GSP coincides with Y for 
d = 2. For every ci > 1, the standard superadditivity argument yields the existence of A^ > 
such that 

F[ud{s,x) < 0, Vs e [l,t]] = t-^'^+^W 

and one has 6 = A2. In [22], an empirical approach using independent interval approximation 
is described, proposing A^ as the first zero on the negative axis of the function 



X ^ 1 + J'^(tix-2x^ f e-"*sin"i[l/(cosh(t/2))^/2]^tJ 



which yields the numerical values Ai = 0.1207, A2 = 0.1862 and A3 = 0.2358. In the above 
formula, the function d 1— )■ A^ is increasing, which is somehow in heuristic accordance with 
the fact that in the first-order expansion 

l/(cosh(t/2))'^/2 = 1 - — + o{t'), 

lb 

the coefficient d 1— t- d/16 also increases. In ^48j it is argued that A^ ~ cyrf at infinity, for 
some constant c > 0. The paper [18] also establishes for every d > 1 a general connection 
between the survival analysis of the equation (14. 5 p and the positivity of a family of random 
polynomials defined as 

n-l 

i=l 

where {^j, i > 0} is a i.i.d. sequence of A/'(0, 1) random variables and X is the deterministic 
variable. Setting A^^ for the number of its real roots, it is argued in [H] that 

(4.6) ¥[Nt+i=Q] = n-2(^2+^'*)+°(^\ 

The method relies on an analysis of the mean density p'^ of the real roots, which is defined 

as 



r-fe 
Td\ 



nN>M] = / Pn{x)dx 
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for every a < b, where N^[a,b] is the number of real roots of P^{X) on [a,b]. It is proved 
in [IB] that the hmit mean density pj^ has a shape which is independent of d on (—1, 1)'^ 
(whence the contribution 2A2 in the exponent), and depends on d on (—1,1) (whence the 
contribution 2Xd). The rate of explosion at ±1 is also explicitly evaluated, which makes 
it possible to obtain the precise link between A^^(— 1, 1) and the process Ud, resp. between 
A''^(— 1, 1)^ and the process M2- We confess not having checked the whole argument of ^8] 
in detail, which seems to deserve a more complete mathematical explanation as it is done in 
[17] for the case d = 2. 

4.3. Wetting models with Laplacian interactions. Let / be a bounded and every- 
where positive probability density over M, centered and having finite variance. Introduce the 
Hamiltonian 'H[a^b]{f), defined for a,b & Z with b — a > 2 and for (f : {a, . . . , 6} — )■ M by 

6-1 



y-laMiv) = Yl ^i^Vn) 



n=a+l 

where V = — log(/) is the potential and 

is the discrete Laplacian on Z. The free pinning model with Laplacian interaction is the 
probability measure on M^~^ defined by 



T^V ( A A \ exp(-H[_i,7v+i](<^)) , , 



where ^qtv i^ ^^^ normalization constant which is called the partition function, and where 
the boundary conditions are given by (/9„i = y^o = ^n = '^N+i = 0. This probability measure 
modelises a certain (1 + 1)— dimensional field (viz. a linear chain {{n,ipn), n = . . . A^}) 
with zero boundary conditions and whose interacting structure is described by the discrete 
Laplacian and the potential V. This chain can be viewed as an example of a discrete random 
polymer in (1 + 1)— dimension. 

The free pinning model with gradient interaction, where A is replaced by the discrete 
gradient Vipn = ^Pn+i — V'n and where the boundary conditions are (po = ipjy = 0, has been 
well studied in the literature and has a natural interpretation in terms of random bridges 
with increment density given by /. The model with Laplacian interaction has exactly the 
same interpretation in terms of integrated random bridges. Specifically, one can easily show 
- see Section 2 in [14] for details - that Pq jy is, with the notations of Section 13. H the 
law of an integrated random walk {A^ = Si + ■ ■ ■ + Sn, n = 1 ... A^ — 1} conditioned on 
An = An^i = 0, with increment density given by /. The partition function Zq ^ is then the 
value at (0,0) of the density of {Sn+i,An+i)- Notice that both above free pinning models 
have natural counterparts in continuous time in the context of semiflexible polymers. The 
gradient interacting case corresponds to directed polymers, whereas the Laplacian interacting 
case corresponds to polymers with non-zero bending energy - see [13] for details. 

The connection with persistence of integrated random bridges is made in considering the 
corresponding wetting model with Laplacian interactions, which is the probability measure 
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on M^~^ defined by 

P^^(dv9i,..., dv9Af_i) = Fl^^{dipi,...,dipN^i\^Pi>0,...,ipN-i>0) 

exp(-?^[_i,iv+i]M), , , 

= ^ l{¥'i>o,...,Viv-i>o} dv?i . . . d(pN-i 

where Z^j^ is the normalization constant, and with the same boundary conditions y^-i = 
<^o = 'Pn = y^N+i = 0. In this model, the discrete random polymer is in the presence of a one- 
dimensional hard wall at zero which forces it to stay non negative. From the very definition, 
one sees that Pjf'^ is the law of the above integrated random walk {An, n = 1 ... A^ — 1} 
whose increment has the density /, and conditioned on fi~}^_i n {A^ = A^^i = 0} with the 
notation ^~}^_i = {Ai > 0, . . . , A^^i > 0}. It is then easily shown that the partition function 
is given by 

Zo^^v = n^N-i I An = Ajv+i = 0]/jv(0, 0) 

where /Ar(0,0) is the value at (0,0) of the density of {Sn+i,An^i). As a consequence of a 
local limit theorem - see Section 2 in [H] for details - it can be shown that /Ar(0, 0) ~ cN'"^ 
at infinity for some explicit constant c > 0. Hence the behaviour of Z^j^ for A^ large, which 
has some importance in physics, is specified by the persistence probability 

F[n+_,\AN = AN+i = 0]. 

The latter quantity has also some independent interest as a question about entropic repulsion 
- see all the references listed in [Hj for more on this subject, and the following is stated in 

m 



Conjecture 13 (Caravenna-Deuschel). With the above notations, one has 

P[fi+_i I An = An+1 = 0] X Ar-V2 

for every centered increment law fi having finite variance. 

This conjecture is related with Conjecture 1 since the event {f^^_i | An = An+i = 0} can 
be decomposed into {Ai > 0, . . . , An/2 > | An = An+i = 0} fl {An/2+i > 0, . . . , An-i > 
I An = An+1 = 0}, the intersection of two roughly independent events with roughly the 
same probability W[Ai > 0, . . . , An/2 > 0], a quantity which should behave like N~^^'^. Notice 
that in the context of semifiexible polymers, a continuous couterpart of ^n-i i^ ^^^ ^^^^ 
when fl is Gaussian was investigated (without conditioning) in [13], where the estimate (13.21) 
is proved. In [14|, the following weak bounds are obtained 

< T>[Q+_^ I An = An+1 = 0] < 



c+ 



AT- - -L-W-l,^- .^.v^x .J _ ^j^g^^ 

for some constants c, C,C- > and c+ > 1. The lower bound entails that the free energy 
vanishes: 

hm ^logZo-^ = 0, 
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whereas the fact that c+ > 1 in the upper bound is crucial to show that the phase transition 
of the wetting model with reward, which is simply the value of the positive parameter e after 
which the free energy of the probability measure 

Af-l 



— ^ J_J_(e(5o(a¥?„) + l{^„>o}d(/?, 



^e,N „=i 



becomes positive, is of first order. In |T3], to which we again refer for more details, it is 
mentioned that Conjecture 13 would yield some further path results for the wetting model 
with reward at criticality. 

4.4. Other physical applications. 

4.4.1. Spatial persistence for fluctuating interfaces. A fluctuating interface is a function h : 
M+ X M'' — )• M evolving in time, with dynamics governed by a certain random equation. The 
problem of spatial persistence concerns the probability p{l) that such a fluctuating interface 
stays above its initial value over a large distance / from a given point in space. One expects 
a behaviour like p{l) = l^^+°W for a positive number 9 independent of the direction, which 
is called the spatial persistence of the interface. In [ID], this question is adressed for the 
Gaussian interface h{t,x) solution to the equation 

(4.7) ^ = -i-Ay/^h + e 

where A is the d— dimensional Laplacian, C, a space-time white noise with zero mean, and 
z > d some fractional parameter, and it is shown with heuristic arguments based on Fourier 
inversion that if p = {z — d + l)/2, then the fractional derivative in any direction xi 

d^h d f I n ,, ,, ,i_„ , 

h{t,y,X2,...,Xd){xi-y) dy 



dx{ 54^1+1 Vr(2- a 

where we have decomposed p = [p] + a into integer and fractional parts, is a one- dimensional 
white noise. This relates the spatial persistence probability p{l) to the persistence probability 
of the Riemann-Liouville process 



1 /■* 
A", = —— / it-s)P-^dBs 

* r(p)yo^ ^ 



In [40] two different regimes are considered. The coarsening one, where the reference point 
is fixed, yields a spatial persistence exponent 6 = 6{p) with the notations of Paragraph 13.31 
The stationary one, where the reference point is sampled uniformly from the ensemble of 
steady state configurations, yields from (13. 5p a spatial persistence exponent 6* = (1 — p)+. 
In the coarsening regime, this entails that the zero crossings of Gaussian interfaces governed 
by (14.71) undergo a morphological transition at z = d + 2, because then p = 3/2. 
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4.4.2. Clustering of sticky particles at critical time. In this last paragraph we consider a 
random walk {Sn, n > 1} with positive increments having expectation E[S'i] = 1. If ra 
particles are fixed at the respective positions i~^Si,i = 1 . . .n with zero initial speed and 
then move according to the laws of gravitational attraction, these particles end up in sticking 
together with conservation of mass and momentum, forming new particles called clusters. 
One is then interested in the number of clusters Kn{t) G [1, . . . , n] viz. the total number of 
particles present at time t > 0. This is a so-called sticky particle model, which is for n large 
connected to the inviscid Burgers equation with random initial data (coupled with some 
scalar transport equation, see [12] for details). This is also an aggregation model having 
connections with astrophysics, and we refer to the introduction of [M] for a clarification of 
these relations and the complete dissipation of all possible misunderstanding. 

The normalization E[S'i] = 1 entails that T — )• 1 in probability, where T stands for the 
random terminal time where all particles have aggregated in a single cluster. A more precise 
result is obtained in [51] in the case when 5*1 has uniform or standard Poissonian distribution, 
namely that the random function 

Knit) - n(l - t^) 



'n 

converges in law to some Gaussian process on the Skorokhod space P[0, 1 — s] for any e > 0. 
In particular, the above quantity converges to some Gaussian law at each fixed time t < 1. 
The situation is however different at the critical time t = 1, at least when 5*1 has a standard 
Poissonian distribution. In this case it can be proved that KniX)/ \/n does not converge to 
zero (the only non-negative Gaussian distribution) as could be expected, and this fact is 
actually a consequence - see [5l] for details - of the estimate 



P 



min y^{Tj-i) > 



n-'l^ 



_ _0 
i=i 

where {r„, n > 1} is a random walk with exponential increments - this latter estimate 
follows from the main result of [I6| . 
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